Boundary critical phenomena of the random transverse Ising model in D > 2 
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Using the strong disorder renormalization group method we study numerically the critical behavior 
of the random transverse Ising model at a free surface, at a corner and at an edge in D = 2, 3 and 
4-dimensional lattices. The surface magnetization exponents are found to be: x B = 1.60(2), 2.65(15) 
and 3.7(1) in D = 2, 3 and 4, respectively, which do not depend on the form of disorder. We have 
also studied critical magnetization profiles in slab, pyramid and wedge geometries with fixed-free 
boundary conditions and analyzed their scaling behavior. 



I. INTRODUCTION 

The quantum Ising model with random couplings 
and/or with random transverse fields (RTIM) is the pro- 
totype of disordered quantum magnets having discrete 
symmetry. This model has a zero-temperature quantum 
phase transition, the properties of which have been stud- 
ied by a special strong disorder renormalization group 
(SDRG) method^. In this method the strongest local 
terms of the Hamiltonian are successively eliminated and 
at the same time new terms are generated perturbatively 
between remaining degrees of freedom^. In one dimen- 
sion (ID) where the topology of the lattice stays invariant 
under the transformation the SDRG equations have been 
solved analytically in the vicinity of the quantum critical 
point 3 -. In this case the phase transition is shown to be 
controlled by a so called infinite disorder fixed point^ 
(IDFP), in which disorder fluctuations are completely 
dominant over quantum fluctuations and therefore the 
renormalization steps are asymptotically exact for large 
scales. Indeed the SDRG results in ID are consistent with 
findings of other analytical^ and numerical methods^—. 

In higher dimensional lattices the topology of the 
lattice is changed during the SDRG steps, therefore 
the SDRG method has to be implemented numerically. 
The first numerical calculations have been performed 
in 2Di£~— and more recently an efficient numerical 
algorithm^— of the present authors made possible to 
extend the calculations^— to 3D and 4D, as well as to 
Erdos-Rcnyi random graphs, which are infinite dimen- 
sional lattices. In all dimensions the phase transition is 
found to be controlled by an IDFP, which justifies that 
the SDRG method provides asymptotically exact results 
for large systems in higher dimensions, too. Quantum 
Monte Carlo simulations for the 2D RTIM are consistent 
with the SDRG results^. Similarly, simulation results 
for the random contact proces o 22 ! 23 - which is expected 
to be in the same universality class^ as the RTIM - arc 
in agreement with the SDRG results in 2D and in 3D. 

In D > 1 dimensions during the SDRG iterations a 
large number of new couplings are generated between re- 
mote sites, which makes the numerical implementation of 
the method rather cumbersome. To avoid this problem 



another, more simple approximation methods have been 
developed and applied to the RTIM^r^i. One of those 2 ^ 
is based on the quantum cavity approach 3 ^, which is 
found to reproduce some of the exact results in ID. How- 
ever, in the Bethe lattice with an effective dimensional- 
ity of D c ff = 2 the method has predicted conventional 
random critical behavior instead of IDFP scaling. The 
quantum cavity method is shown to be equivalent to a 
linearized transfer matrix approach^!. If no linearization 
is performed (this is the so called non-linear transfer ap- 
proach) than the method has lead to IDFP behavior for 
D > 2, too2i. Also approximate renormalization group 
schemes have been suggested^— , during which the or- 
der of the RG steps is changed in such a way that the 
proliferation of new couplings is avoided. These methods 
have reproduced some exact ID results and also provide 
IDFP behavior for D > 2, in agreement with the stan- 
dard SDRG method. 

Most of the results about the critical behavior of the 
RTIM have been calculated for bulk quantities. For ex- 
ample the order-parameter of the RTIM is the average 
magnetization and its value in the bulk, mj, has the scal- 
ing behavior m\> ~ L~ x , where L is the linear size of the 
system and x is the scaling exponent of the bulk magne- 
tization. Real systems, however, have finite extent and 
they are limited by boundaries. At a free surface the 
scaling behavior of the average surface magnetization, 



m s , involves a new exponen' 
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Due to missing 



bonds at the surface there is weaker order, therefore gen- 
erally x s > x. For the ID RTIM several properties of 
the surface magnetization (the distribution function, av- 
erage and typical behavior, etc) is exactly known^^. 
For example the surface scaling exponent, x s = 1/2, is re- 
lated to the persistence properties of ID random walks^. 

For higher dimensional RTIM less attention is paid 
to the calculation of the surface magnetization: we are 
aware of one recent work 2 ^, in which the surface magne- 
tization exponent has been calculated by the non-linear 
transfer approach. The obtained values are x s = 1.2 and 
1.34, in 2D and 3D, respectively, which are to be com- 
pared with the SDRG results for the bulk magnetization 
exponent^—: x = 0.98 and 1.84, in 2D and 3D, respec- 
tively. Since in 3D the surface magnetization exponent 
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of the non- linear transfer approach is smaller, than the 
expected correct value of the bulk magnetization expo- 
nent we conclude that the non-linear transfer approach 
underestimates the values of x s . Therefore there is a ne- 
cessity to obtain more accurate estimates for the surface 
critical properties of the RTIM. 

In this paper we study the boundary critical behavior 
of the RTIM in higher dimensional systems for D = 2, 3 
and 4 by the SDRG method. In the calculation we use 
the numerical algorithm, which has been developed in 
Refsi 17 ' 18 and has been used to study the bulk critical 
behavior of the systems. Besides the surface magneti- 
zation exponent, x a , we calculate local magnetization 
exponents^ 7 -, which are associated with corners (in 2D 
and 3D) as well as with edges (in 3D). We also calcu- 
late critical magnetization profiles, when spins are fixed 
at some surfaces of the system and study their scaling 
properties. 

The structure of the paper is the following. The 
essence of the SDRG method and its application to the 
calculation of the boundary magnetization is described 
in SecHU Our results for 2D, 3D and 4D lattices are pre- 
sented in Scc lIIII and discussed in the final section. 



II. SDRG CALCULATION OF THE LOCAL 
MAGNETIZATION 

Here we consider the boundary critical properties of 
the RTIM defined by the Hamiltonian: 



H 



(i) 



in terms of the a^' z Pauli operators at site i of a D- 
dimensional cubic lattice. The Jij > nearest-neighbor 
couplings and the hi > transverse fields are indepen- 
dent random numbers taken from the distributions p( J) 
and q(h), respectively. In this paper we use two different 
disorder distributions, which have already been used in 
Refs for the calculation of the bulk properties. The 
advantage of using these distributions is, that the loca- 
tion of the critical points have already been determined. 
In the bulk calculations we have observed that the two 
distributions lead to identical critical exponents, within 
the error of the numerical method. Here we assume that 
universality holds to local critical properties, too, and 
check this assumption numerically. In both type of dis- 
order the couplings arc taken from a uniform distribution: 
p(J) = 0(J)0(1 — J), where Q(x) is the Hcaviside step- 
function. For the box-h model the transverse fields have 
a box-like distribution q(h) = ^-0(ft,)0(/i& — h), whereas 
for the fixed-h model the transverse fields are constant: 
q{h)=S(h-hf). 

In the calculations of different local magnetizations 
(surface, edge and corner) we have used different finite ge- 
ometries, in which fixed spin boundary conditions (b.c.) 
have been used at given planes and the magnetization 
profile, mi, is measured perpendicular to the fixed planes. 



(These are located at I = 1, thus mi = 1). For the sur- 
face magnetization a slab of size L x N ^ 1 (L < N) is 
used, and in the short direction we use fixed- free b.c, 
while in the other (D — Indirections periodic b.c. is ap- 
plied. The surface magnetization is given by m s = mi = L, 
which scales at the critical point as m s ~ L~ Xs . 

The corner magnetization is measured in 2D at the free 
corner of a half square and in 3D at the free corner of a 
cube, having the shape of a pyramid. (In the following we 
use the term pyramid in the 2D case, too.) The spins at 
the base of the pyramid are fixed, while at other surfaces 
free b.c.-s are used. The magnetization profile, mi is 
measured between the base (I = 1) and the corner (I = 

L c = ^p-L), and the corner magnetization is given by 
m c = mi=L e - This scales at the critical point as m c ~ 
L~ Xc with the corner exponent, x c . 

Edge magnetization is calculated in 3D at the free edge 
of length N of a square column of size LxLxN,L<N. 
The square column is cut at the square diagonal plane, 
thus have the shape of a wedge, and the spins at the base 
of the wedge are fixed. In the long direction periodic 
b.c. is used, while at the other two symmetric surfaces 
free b.c. is applied. The magnetization profile, m; is 
measured between the base of the wedge (1 = 1) and the 

free edge (I = L e = ^rL). mi is translationally invariant 
along the long direction and the edge magnetization is 
given by m e = mi—L e - This scales at the critical point 
as m e ~ L~ x ", with the edge exponent, x e . The applied 
geometries in 3D are shown in Fig. [TJ 
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FIG. 1: (Color online) The three geometries used in the 
calculation of the local magnetization in 3D. Left panel: slab; 
middle panel: wedge; right panel: pyramid. Spins at the 
shaded planes are fixed and the arrows indicate the directions 
in which the density profiles are measured. 

We note that the corners and edges we consider here 
have the specific opening angle: 7r/2. The local critical 
exponents are generally angle dependent^ 7 -, but we do not 
study this problem in the present paper. 

To calculate the local magnetization in the different 
geometries we have used the SDRG method, which is an 
iterative procedure working in the energy space. At each 
step the largest local term of the Hamiltonian, either a 
coupling, Jij, or a transverse field, hi, is decimated and 
new terms are generated between the remaining sites in 
a perturbation calculation. For coupling decimation the 
two sites, i and j with original magnetic moments, /ii and 
Hj , are merged to a new cluster with an effective moment 
/4j = l l i + Mjj which is placed in an effective transverse 
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field of strength: = hihj/Jij. In transverse field dec- 
imation the site i is eliminated and its nearest-neighbor 
sites, say j and k will be connected by an effective cou- 
pling: Jj k = ma,x{JjiJik/hi, Jjk}- In this last step the 
so called maximum rule is applied, the use of which is 
justified at an IDFP. 

We apply the numerical algorithm of the SDRG 
method in Rcfsi 17 ' 18 , which has been used to locate the 
critical point of the system (for the two forms of the dis- 
order) and to calculate the bulk critical exponents at the 
IDFP for different dimensions, D = 2, 3 and 4 (which are 
disorder independent and listed in Table Hi]). To calculate 
the critical magnetization profile, mi , we renormalize the 
system up to the last effective site and consider the effec- 
tive cluster, C, which contains the fixed sites at I = 1. If 
the system has an IDFP, then all spins of C are strongly 
correlated: in leading order all these spins point to the 
same direction as at I = 1, whereas other sites (not con- 
tained in C) have negligible contribution to the longitu- 
dinal magnetization. Let us denote by ni the number of 
sites in C at position I and the number of equivalent sites 
by n (it is N d ~ x , N and 1 in the slab, wedge and pyramid 
geometry, respectively). The average value of the local 
magnetization is then given by: mi — [ni/h] av , where 
[. . . ]av stands for the average over disorder realizations. 

At the critical point the asymptotic form of the mag- 
netization profile is given by scaling considerations. Ac- 
cording to Fisher and de Gennes^ the decay of the mag- 
netization from the fixed surface is given by: 

mi~l~ Xb , x = x b , KZ<L, (2) 

thus it includes the bulk magnetization exponent. Close 
to the free endpoint (surface, corner or edge) the magne- 
tization profile has a different power-law decay^: 

m r ~ (Z') Xq6 , l«r = I a -; + l«L a , (3) 

with x a b = x a — x and a relates to the type of endpoint : 
s, c or e (and L s = L). These relations will be used to 
obtain independent estimates for the local magnetization 
exponents. The two scaling relations in Eqs.@ and (O 
can be incorporated into an interpolation formula: 

mi = A [ sm (7rA)r [cos(7rA/2)]*" , (4) 

with A = This relation is exact for ID conformally 
invariant quantum systems^. Although the RTIM is not 
conformally invariant, in ID Eq.@ is found to be an 
excellent approximation^ .. In the following in the slab 
geometry (a = s) we shall check the accuracy of Eq.(j4]) 
in higher dimensions, too. 

III. RESULTS 

We have calculated the magnetization profiles in the 
three geometries described in Scc|TI] in different dimen- 
sions: 2 < D < 4 by the SDRG method using two differ- 
ent forms of disorder. The largest sizes of the systems, 



the typical aspect ratios of slabs and wedges as well as 
the typical number of disorder realizations are collected 
in Table |H Since only a small fraction of samples con- 
tains such a correlation cluster, C, which have also sites 
at the free extremity of the system (surface, edge or cor- 
ner) one should consider a large number of realizations. 
For surfaces and edges in a given sample there are several 
end-point positions, for which we perform the averages. 
For corners, however, there is just one end-point in a 
sample, therefore one should take even larger number of 
realizations. In the following we present our numerical 
results obtained in different dimensions. 

TABLE I: Details of the numerical calculation of the local 
magnetization. L max : largest linear size; N/L: typical aspect 
ratio; iV # : typical number of realizations. 
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A. Calculations in 2D 

1. Surface magnetization 
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FIG. 2: (Color online) Magnetization profiles in 2D in 
the slab geometry for fixed-free b.c.-s in a system of width 
L = 512 for box-h and fixed-/i randomness. The interpola- 
tion formula in Eq.Q is represented by dashed lines. In the 
inset the ratio of the magnetization profile and the interpo- 
lation formula in Eq.Q is shown for x = 0.982, x s — 1.6, 
A flxcd = 1.20 and A hox = 0.282. 

The magnetization profiles in the slab geometry calcu- 
lated by the two types of disorder are shown in FigJ5]as 
a function of the relative position: A = l/L, see Eq.@. 
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Here we use a finite-size shift of lo = 0(1) at the bound- 
aries. As already observed in the calculation of the bulk 
magnetization the typical correlation clusters for fixed- h 
disorder contain approximately 6-times more sites, than 
for bax-h disorder. As a consequence the magnetization 
profiles are also comparatively larger for fixed-/i disor- 
der. As seen in Fig 12] the magnetization is monotonously 
decreasing and the variation is very fast near the two 
end-points, which are then analyzed in log-log plots in 
Figsj3]and|H respectively. 
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FIG. 4: (Color online) Magnetization profiles near the free 
boundary in 2D for the slab geometry for the two type of 
randomness with L — 512. In both cases the decay is char- 
acterized by the exponent, x s t, — 0.65(2). In the inset the 
finite-size estimates for the surface magnetization exponent 
are presented. The extrapolated (disorder independent) value 
is given in Table [TT1 



FIG. 3: (Color online) Magnetization profiles near the fixed 
boundary in 2D for the slab and the pyramid geometries for 
the fixes-/i and box-h randomness with L = 512 and 256, 
respectively. In all cases the decay is characterized by the 
same exponent, Xb = 0.98(1), which according to the Fisher- 
de Gennes result in Eq.([2| is equivalent to the bulk magneti- 
zation exponent, see Table ITT1 



formula, in which the exponents in Table |TT] have been 
used. As seen in this figure the interpolation formula 
represents a good approximation, but the agreement is 
not perfect, the largest discrepancy is about 10%. 



Close to the fixed boundary the magnetization pro- 
files are shown in Fig|3l together with the similar pro- 
files in the pyramid geometry, which will be analyzed 
in Sec IIII A"2] The magnetization profiles for the two 
disorder have the same power-law decay and the de- 
cay exponent is estimated from the largest systems as 
= 0.98(1). This is to be compared with the value of 
the bulk magnetization exponent x = 0.982(15), which 
has been calculated in RefJ£ by finite size scaling. We 
can thus conclude that the Fisher-de Gennes scaling pre- 
diction in Eq.@ is well satisfied. 

Also at the free-boundary the profiles have a power-law 
variation (see Figj4j) and the corresponding exponent is 
estimated from the largest system as: x sb — 0.65(2). We 
have also calculated the surface magnetization exponent, 
x Sl from the finite-size scaling behavior of the surface 
magnetization, m s . Two-point estimates for x s are pre- 
sented in the inset of Fig. 2J which have the same limiting 
value for large L for the two type of disorder, which is 
presented in Table [TTJ Comparing x sb with x s — x we can 
conclude that the scaling prediction in Eq.Q is satisfied. 

We have also checked the accuracy of the interpolation 
formula in Eq. (0| and in the inset of FigJ5]we have plotted 
the ratio of the measured profile and the interpolation 



2. Corner magnetization 

The calculations are performed in the pyramid geom- 
etry and the critical magnetization profile close to the 
fixed plane is shown in Fig[2] in a log-log plot for the 
two different initial disorder. As discussed in Scc lHI ATI 
in this figure also the profiles in the slab geometry are 
presented and the two types of profiles are very close to 
each other: they are indistinguishable within the error of 
the calculation. Thus in agreement with scaling theory 
the decay of the profile in the pyramid geometry is in 
a power-law form with a decay exponent, x b — x. The 
magnetization profile at the other end, i.e. starting from 
the corner is shown in Fig J5] and the corresponding decay 
exponent of the magnetization, x c b, is presented in Table 

m 

From finite-size scaling the corner magnetization expo- 
nent, x c , is calculated by two-point fit and the effective, 
size-dependent exponents are presented in the inset of 
Fig JS] for the two different type of disorder. The extrap- 
olated value which is disorder independent is given in 
Table IE 
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FIG. 5: (Color online) Magnetization profiles in the pyramid 
geometry near the corner in 2D for the two type of randomness 
with L = 256. In both cases the decay is characterized by the 
exponent, x c t = 1.35(10). In the inset the finite-size estimates 
for the corner magnetization exponent are presented^. The 
extrapolated (disorder independent) value is given in Table 



B. Calculations in 3D 

1. Surface magnetization 



exponent, see Tabic |TTJ Near the free surface the profiles 
for the two types of disorder are shown in Fig. [JJ and the 
estimated decay exponent, x s b, is presented in Table [TT1 
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FIG. 7: (Color online) Magnetization profiles near the free 
boundary in 3D for the slab geometry for the two types of 
randomness with L = 64. In both cases the decay is char- 
acterized by the exponent, x s b = 0.84(7). In the inset the 
finite-size estimates for the surface magnetization exponent 
are presented. The extrapolated (disorder independent) value 
is given in Table HTl 



The surface magnetization exponent is estimated 




through finite-size scaling and the effective, size- 
dependent values are shown in the inset of Fig. [7] for 
the two types of disorder. The extrapolated exponent is 
disorder independent and given in Table [TTJ We conclude 
that the scaling relation in Eq.Q is satisfied within the 
error of the calculation. 

Wc have checked the accuracy of the interpolation for- 
mula in Eq.Q and the ratio of the measured profile and 
the interpolation formula is shown in the inset of Fig. 
|6l Also in this case Eq.Q is a good approximation, the 
maximal discrepancy is somewhat larger, than in the 2D 
case, see in Fig. [5J 



FIG. 6: (Color online) Magnetization profiles in 3D in the 
slab geometry for fixed-free b.c.-s in a system of width L = 64 
for box-h and fixed-/i randomness. The interpolation formula 
in Eq.(|4]) is represented by dashed lines. In the inset the ratio 
of the magnetization profile and the interpolation formula in 
Eq.(g]) is shown for x = 1.84, x s = 2.65, ^lfl xe d = 1-78 and 
A box = 0.316. 

The magnetization profiles in the slab geometry are 
shown in Fig. [6] for the two types of disorder. Close to 
the fixed boundary the exponent associated to the decay 
of the profile is estimated as xi = 1.855(20) which agrees 
with the finite-size estimate of the bulk magnetization 



2. Edge magnetization 

We have measured the magnetization profile in the 
wedge geometry and here we analyze its behavior close 
to the free edge, see Fig. [H 

The estimated decay exponent, x e b, is presented in Ta- 
blc|TT]together with the extrapolated value of the edge ex- 
ponent, x e , for which the finite-size estimates are shown 
in the inset of Fig. [5] In this case, too the scaling relation 
in Eq.Q is satisfied. 
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FIG. 8: (Color online) Magnetization profiles near the free 
edge in 3D in the wedge geometry for the two type of random- 
ness with L = 64. In both cases the decay is characterized by 
the exponent, x e b — 1.75(15). In the inset the finite-size es- 
timates for the edge magnetization exponent are presented^. 
The extrapolated (disorder independent) value is given in Ta- 
ble HH 
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FIG. 9: (Color online) Magnetization profiles near the free 
corner in 3D in the pyramid geometry for the two type of ran- 
domness with L = 64. In both cases the decay is characterized 
by the exponent, x c b = 2.65(25). In the inset the finite-size 
estimates for the edge magnetization exponent are presented 
up to L = 32^i. The extrapolated (disorder independent) 
value is given in Table ITTl 



3. Corner magnetization 

We close our study in 3D by calculating the magnetiza- 
tion profile in the pyramid geometry: the result is shown 
in Fig. IHl close to the free corner. (In the inset finite-size 
estimates of the corner exponent are presented.) Esti- 
mates of the decay exponent, x c b, and the corner ex- 
ponent, x c , are presented in Table ITTl which satisfy the 
scaling relation in Eq.Q. 



C. Calculations in 4D 

In 4D the available system sizes are limited, see Table 
UJ therefore we could only study the magnetization profile 
in the slab geometry, which is shown in Fig. [TO] 
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FIG. 10: (Color online) Magnetization profiles in 4D in the 
slab geometry for fixed-free b.c.-s in a system of width L = 32 
for box-h and fixed-/i randomness. The interpolation formula 
in Eq.Q is represented by dashed lines. In the inset the ratio 
of the magnetization profile and the interpolation formula is 
shown for x — 2.72, x 3 = 3.7, Afl xc d = 4.19 and A^ ox = 0.625 



Close to the fixed surface the decay exponent is cal- 
culated as Xb = 2.72(10), which agrees well with the 
finite-size estimate of the bulk magnetization exponent, 
see Table |TT] The magnetization profile close to the free 
surface is shown in Fig. 1111 

Estimates of the decay exponent, x s b, and that of the 
surface magnetization exponent, x Sl which are presented 
in Table ITTl contain somewhat larger errors, than in lower 
dimensional calculations. However the scaling relation in 
Eq.fJH]) is satisfied in this case, too. Also the interpolation 
formula in Eq. @ is a good approximation as can be seen 
in the inset of Fig. [TO] 



IV. DISCUSSION 

In this paper we have used the SDRG method to calcu- 
late the magnetization profiles of the random transverse 
Ising model in 2D, 3D and 4D in different geometries: 
slab, corner and wedge having a fixed surface. At the 
critical point decay exponents are calculated both at the 
fixed end and at the free end of the profiles. These ex- 
ponents, which are presented in Table UTI are found to be 
disorder independent, at least for strong enough disor- 
der, for which the critical properties of the system are 
controlled by an IDFP. From finite-size scaling studies 
of the local magnetization at the free ends of the profile 
local (surface, corner and edge) critical exponents are cal- 
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FIG. 11: (Color online) Magnetization profiles near the free 
boundary in 4D for the slab geometry for the two type of 
randomness with L = 32. In both cases the decay is char- 
acterized by the exponent, x s b = 0.85(15). In the inset the 
finite-size estimates for the surface magnetization exponent 
are presented. The extrapolated (disorder independent) value 
is given in Table [Til 

culatcd, see Table [TT1 For all types of profiles considered 
here the scaling relations in Eqs.® and © arc satisfied 
and the interpolation formula in Eq.(j4|) is found to be a 
good approximation in the slab geometry. Our results, 
concerning the properties of the average local magnetiza- 
tion of the RTIM are rather complete, these are compa- 
rable with the existing results in the non-random system. 

By the SDRG method the average local magnetiza- 
tion, mi, is obtained as the ratio of such rare realiza- 
tions, in which the correlation cluster contains the given 
site. By this method the typical value of the local mag- 
netization could be estimated by the strength of the 
effective coupling, J' u , which is generated between the 
fixed surface and the site. For surface spins it scales as 
m t s vp ~ J' 1L ~ exp(-AL^), where ip is a characteris- 



tic exponent in the IDPF— ii, which has been calculated 
in RefsJ£~— . Thus by calculating m' yp by some other 
means one can obtain independent estimates for the ex- 
ponent tp. 

Concerning the dimensional dependence of the average 
surface magnetization exponent, we write it in the form: 
x s (D) = D s + p(D), where D s = D — 1 is the dimen- 
sion of the surface and p(D) is a number close to 1/2. In 
ID p(l) — 1/2 is shown to be the persistence exponent 
of the random walk^ and we propose here an analogous 
explanation for D > 1, too. Let us denote by fj, s the num- 
ber of surface points of the correlation cluster, C which 
starts at the fixed boundary. We have checked, that /z s 
has an exponential distribution: -P(/i s ) ~ exp(— fj, s /fi), 
with fi w B{N/L) D - 1 , c.f. in 2D we have B = 0.5 and 
B = 1.0, for Rxcd-h and hox-h disorder, respectively 
Consequently the surface points of C arc grown from 
uncorrclated domains. The average number of surface 
points in an area L u ~ x then scales as [/x s ] a v ~ L~ p ( D \ 
As explained in RefJ^ (see Fig. 5 there) the correlation 
cluster is embedded into a connected subgraph, which 
contains all the decimated sites (the results of both h- 
and J-dccimations) and which is related to a low-energy 
excitation of the system. The number of points in the 
connected subgraph is L ~ L Df , where the fractal di- 
mension, Df, is close to one. If we replace the connected 
subgraph with a linear chain with L sites then we obtain 
from the random walk result: p(D) « -D//2. Indeed our 
numerical estimates^ of Df and the surface magnetiza- 
tion scaling dimensions in Table [TT] are in agreement with 
this relation. In higher dimensions we expect that the 
structure of the SDRG transformation, in particular the 
topology of the connected clusters follows the trend ob- 
served in this paper, thus the surface magnetization expo- 
nent generally obeys the relation: x s (D + l) —x s {D) « 1. 
Based on this result we expect that a simplified SDRG 
procedure, such as a modified version of those used in 
Refs3£~— can be constructed, which captures the main 
results about the surface critical properties of the RTIM. 



TABLE II: Estimates of the critical exponents obtained by 
finite-size scaling: x, x„, x c , x e and the exponents associated 
with the decay of the profile: Xb, x s t, x c b, x e b- (x is taken from 
Refs.— ~— and the exact results in ID are fromRefs.-^^). For 
the pure system (see c.f. ir*^) x s = 0.5 and 1.27 for ID and 
2D, respectively and for D > 3 the mean-field result holds: 
x s = (D + l)/2. The corner exponent in 2D is x c — 2.06. 





bulk 


surface 


corner 


edge 




X 


Xb 


X s 


Xsb 


X c 


Xcb 


x e 




ID 


(3 - v/5)/4 




0.5 












2D 


0.982(15) 


0.98(1) 


1.60(2) 


0.65(2) 


2.3(1) 


1.35(10) 






3D 


1.840(15) 


1.855(20) 


2.65(15) 


0.84(7) 


4.2(2) 


2.65(25) 


3.50(15) 


1.75(15) 


4D 


2.72(12) 


2.72(10) 


3.7(1) 


0.85(15) 
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Our results about the local critical behavior of the 
RTIM are relevant to other random quantum magnets 
having discrete symmetry, we mention the random quan- 
tum Potts^, clock and Ashkin- Teller models^. Also the 
surface, corner and/or edge exponents of the random con- 
tact process are expected to be given by the RTIM values 
in Table [TTJ To check this conjecture one should repeat 
recent Monte Carlo simulations about this mode l 22 i 23 . 

Our studies of the local critical behavior can be ex- 
tended in different directions. For example, one can mea- 
sure the corner and edge magnetization exponents at dif- 
ferent opening angles or one can consider anisotropic sys- 
tems, in which the distribution of disorder is different in 



the different directions. One can also consider the sur- 
face critical behavior in the presence of enhanced surface 
couplings, in which case the so called extraordinary and 
surface transitions 3 ^— could be studied, too. 
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